Magnetic excitations in coupled Haldane spin chains near the quantum critical point 
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Two quasi- 1-dimensional S = 1 quantum antiferromagnetic materials, PbNi2V208 and 
SrNi2V2 08, are studied by inelastic neutron scattering on powder samples. While magnetic inter- 
^ . actions in the two systems are found to be very similar, subtle differences in inter-chain interaction 

strengths and magnetic anisotropy are detected. The latter are shown to be responsible for quali- 
tatively different ground state properties: magnetic long-range order in SrNi2V208 and disordered 
"spin liquid" Haldane-gap state in PbNi2V208. 



O 
O 
O 



o 

\ I. INTRODUCTION 

(N 
> 

' After two decades of intensive theoretical and experimental studies, 1-dimensional (ID) Heisenberg quantum an- 
\ tiferromagnets (AF) are now rather well understood. A great deal of work has been done on integeri^gin systems 
CO that have a spin-liquid ground state and a famous Haldane gap in the magnetic excitation spectrum.lila The focus 
in quantum magnetism has now shifted towards studies of more complex phenomena, that include inter-chain inter- 
actions, spin-lattice coupling, and/or spin-vacancies and substitutions. Of particular current interest is the quantum 
phase transition between spin-liquid (non-magnetic) and ordered states. This type of transition in gapped ID sys- 
tems occurs as 3D magnetic interactions and/or magnetic anisotropy are increased beyond certain threshold values. 
^ ■ Their effect is to lower the energy of excitations at certain points in reciprocal space, and ultimately induce to a 
' soft-mode transition|-to,a Neel-like ordered structure. An example of such behavior is found in the extensively studied 
I I CsNiCla compjOund.H^El The corresponding phase diagram has been worked out by several authors, including Sakai 
; and Takahashitl (Fig. |l|). 

, The most direct way to observe such a transition experimentally is by looking at a series of isostructural compounds 
with slightly dilferent inter-chain coupling constants or anisotropy terms. The problem is that most known quasi-lD 
S = 1 AF materials are otherwise deep inside the spin- liquid area of the phase diagram (goQfi.l-D systems), or 
obviously in the 3D Neel-like or XY-like ordered phases (Fig. |lj). CsNiCls and related compoundsOufl are perhaps the 
only systems close to the phase boundary that have been extensively studied to date. Unfortunately, these compounds 
order in 3 dimensions at low temperature, and have no isostructural counterpart with a spin-liquid ground state. Only 
about a year ago the first quasi-l-D integer-spin AF that is still in the spin-liquid state, but is on the verge of a 3D 
ordering instability, was characterized.EI This Haldane-gap material, PbNi2V208 is so dose to the phase boundary 
that LRO, absent in the pure compound, can be induced by spin-vacancy substitution.Ll Moreover, an isostructural 
undoped system, namely SrNi2V208 does order in 3-dimensions at low temperatures. For the first time two very 
similar stoicheometric quasi-lD materials with such vastly different ground state properties can be compared in 
experimental studies. 

The magnetic properties of both PbNi2V208 and SrNi2V208 are due to spin S = 1 octahedrally-coordinated 
Ni^+ ions, while the sites are presumed to be non-magnetic. The crystal structure, visualized in Fig. |^a, 
is tetragonal, space group Mlcd, with lattice constants a = 12.249(3) A, c = 8.354(2) A for PbNi2V208,E3 and 
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a = 12.1617 A, c = 8.1617 A for SrNi2V208," respectively. The magnetic Ni-sites are arranged in peculiar spiral- 
shaped chains that run along the unique crystal axis, as shown in Fig. Even though all nearest neighbor Ni- 
Ni bonds are crystallographically equivalent, the spin-spirals have a step-4 periodicity. The dominant magnetic 
interaction is antiferromagnetic, between nearest-neighbor Ni^+ spins within each chain. The corresponding exchange 
constant, J ~ 8.2 meV in both systems, was deduced from the high-temperature part of the experimental x(T) 
curves.EI The ground state of PbNi2V208 is a Haldane singlet, and the excitation spectrum has an energy gap, as 
unambiguously shown in low-temperature x(T) and C(T) measurements. The energy gaps, 1.2 meV and 2.2 meV, 
for excitations polarized along, and perpendicular to the chain-axis, respectively, were accurately determined in high- 
field magnetization studies. The observed anisotropy of the spin gap is attributed to single-ion easy-axis magnetic 
anisotropy on the Ni-sites D « —0.23 meV. 

Unlike PbNi2V208, SrNi2V20g orders magnetically in three dimensions at Tn = 7 K. The magnetic structure 
has not been determined to date, but, according to bulk measurements, is of a weak-ferromagnetic type, with a 
dominant antiferromagnetic component. The ordered staggered moment is along the unique c crystallographic axis. 
The weak-ferromagnet distortion of this Neel (coUinear) spin arrangement is attributed to the presence of weak 
Dzyaloshinskii-Moriya off-diagonal exchange interactions in the non-centric crystal. 

Preliminary inelastic neutron scattering studiesO provided an estimate for the inter-chain interaction strength. Only 
a limited amount of neutron data are available for PbNi2V208, and, to date, none for SrNi2V208. The present paper 
deals with more extensive comparative inelastic neutron scattering studies of both materials. Our results reveal the 
subtle differences between the two systems, responsible for their vastly distinct ground state properties. 

II. MAGNETIC INTERACTIONS 

Prior to reporting our new experimental findings, we shall briefly discuss the magnetic interactions that may play 
an important role in the physics of PbNi2V208 and SrNi2V208, and construct a model spin Hamiltonian for these 
systems. In doing so, we shall introduce the notation used throughout the rest of the paper. 

As mentioned in the introduction, all previous studies point to that the dominant magnetic interaction is the 
antiferromagnetic coupling J between nearest-neighbor spins (2.78 A) within each spiral-shaped chain. It is also 
clear that taking the expectedly weaker inter-chain interactions into account is crucial to understanding the static and 
dynamic properties. In our previous work (Ref. ^ we assumed that the dominant inter-chain coupling is between pairs 
of Ni^+ ions from adjacent chains, such that both interacting spins have the same c-axis fractional cell coordinate, and 
that each Ni-site is coupled to 4 other sites. The resulting coupling topology is schematized in Fig. Considering 
only these inter-chain bonds, however, can not result in a very realistic model. First of all, they do not correspond to 
the shortest inter-chain Ni-Ni distance. More importantly, in the crystal structure of PbNi2V208, there seem to be no 
obvious superexchange pathways corresponding to these links. If one pays close attention to crystal symmetry, bond 
length, and possible superexchange routes involving the VO4 tetrahedra, one arrives at a more complex interaction 
geometry, visualized in Fig. ^. Here one assumes interactions along the shortest inter-chain Ni-Ni bond (5.0 A). 
Coupled Ni^+ ions are offset relative to each other by c/4 along the chain axis, and are bridged by a VO4 tetrahedron. 
All such bonds are crystallographically equivalent. Each Ni-site is linked with only two adjacent chains. The overall 
coupling topology for this model is as shown in Fig. We shall denote the corresponding exchange constant as Ji. 

Magnetic anisotropy is clearly manifest in bulk susceptibility and magnetization measurements, and must be ex- 
plicitly included in the spin Hamiltonian. Because of the strong dispersion along the chain-axis, all the low-energy 
response of each chain is concentrated at wave vectors close to the 1-D AF zone-center. In this long-wavelength limit, 
2-ion anisotropy of in-chain interactions and single-ion anisotropy associated with individual spins can not be distin- 
guished. In our model we shall therefore include only the latter, and write the corresponding term in the Hamiltonian 
as: 

ffsingle-ion = D'^{S-j,)^. (1) 

i,k 

The choice of the anisotropy axis along the chain direction is based on previous bulk magnetic studies. Here D is the 
anisotropy constant {D < — easy-axis) and Si^k is the spin operator for site i in chain k. At the same time we shall 
assume in-chain exchange interactions to be isotropic and write this Heisenberg term as: 

-f^in— chain = J ^ ^ 'S'l.fe ■S'i-fi^fc . (2) 

i,k 

The magnitude of dispersion perpendicular to the chain axis is expected to be rather small, so the entire range of wave 
vector transfers in the (a, 6) plabe will influence the low-energy properties. For this reason, two-ion anisotropy of Ji, 
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unlike that of J, is a relevant parameter and should be considered. The inter-chain coupling term in the Hamiltonin 
will thus have the form: 



-^intcr— chain 

,k.k 



Here the sum is taken over pairs of next-nearest-neighbor spins. The three terms discussed above constitute the spin 
Hamiltonian that we propose for PbNi2V208 and SrNi2V208. Of course, other magnetic interactions can be active 
in the system as well. As will be discussed below, however, the described model can reproduce the experimental data 
with a minimal number of parameters. 



III. EXPERIMENTAL 



Inelastic neutron scattering studies of PbNi2V208 and SrNi2V208 powder samples (about 10 g each) were carried 
out in two series of experiments. The conventional 3-axis technique was used to measure certain characteristic 
constant- -E and constant-Q scans. These measurements were performed at the TASP and Druchal spectrometers at 
the continuous spallation source SINQ at Paul Scherrer Institute. Neutrons of a fixed final energy of 5 meV or 8 meV 
were used with pyrolytic graphite (PG) monochromator and analyzer, and a typical (open) — 80' — 80' — (open)' 
collimation setup. For some scans a horizontally focusing analyzer was used to increase the useful scattered intensity, 
in which case no collimators were used after the sample. In most of the measurements a flat analyzer was employed 
instead. In the 5 meV-final configuration a Be filter was inserted after the sample to suppress higher-order beam 
contamination. The spectrum of neutrons emerging from the guide is such that no higher-order filter was required 
for the 8 meV-final measurements. 

In the second series of experiments we took advantage of the area-sensitive detector (ASD) 3-axis setup available 
for the NG-5 "SPINS" spectrometer installed at the National Institute of Standards and Technology. In special cases 
when a large domain ol E — Q space needs to be surveyed in a spherically symmetric sample (such as powder) this 
technique can provide an amazing, almost 10-fold increase of data collection rate, as compared to the standard 3-axis 
geometry, without any penalty in resolution. In our experiments on PbNi2V208 and SrNi2V208 we utilized a PG 
monochromator and a (open) — 80' — 80' — 80' (radial) array of collimators. The measurements were done in the fixed- 
final-energy mode, with the central blade of the composite analyzer tuned to — 3.125 meV. In the experiment 
scattering events with final neutron energies in the range Ef ± 0.4 meV are registered simultaneously. To suppress 
higher-order beam contamination we used a Be-0 filter after the sample. The data were taken for momentum and 
energy transfers of up to 2.5 A"-'^ and 7 meV, respectively. 

In both series of experiments the sample environment was a standard "ILL-orange" cryostat, and the temperature 
range 1.5 — 30 K was covered. The background was measured by repeating some scans with the analyzer moved 
away from its elastic position by 10°. In the SPINS experiment the background due to air scattering was separately 
measured with the sample removed from the spectrometer. 



IV. EXPERIMENTAL RESULTS 



A. Results obtained at low temperatures 



The inelastic neutron scattering intensities measured in PbNi2V208 and SrNi2V208 at T = 1.2 K with the ASD 
setup are shown in the false-color plot in Fig. ^ (background subtracted). The resolution of the area-sensitive detector 
is greater than the actual energy and wave vector resolution of the spectrometer (0.11 meV and 0.023 at fiu = 0). 
The customary procedure is to re-bin the mesh data to a coarser grid. Instead, for visualization purposes, the data 
in Fig. 1^ were smeared using a fixed-resolution 0.5 meVx0.06 A^^ FWHM Gaussian filter. As reported previously,H 
constant-energy scans measured in PbNi2V208 at 4 meV energy transfer have a peculiar and very characteristic shape. 
Such scans were collected at T = 2 K for both materials, using the 8 me V- final/flat analyzer setup SINQ. Theses 
scans are shown in Fig. H. In addition, representative constant-Q scans were also measured in this configuration, and 
are shown in Figs. ^ and||. 

Despite the different ground states, the main features of the powder-averaged dynamic cross section for the two 
materials are quite similar. This, however, is not particularly surprising. The dynamic structure factor S{Q,uj) is 
severely smeared out by the spherical averaging that occurs in a powder sample. The main features are defined by 
the form factor of each spiral-shaped chain, and by the steep dispersion along the chain axis, expected to be almost 
identical in the two systems. The significant difference between the two materials is expected to be in the magnitude 
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of the inter-chain interactions, and possibly single-ion anisotropy. These effects are much more difficult to observe, 
as they only influence the weak spin wave dispersion in the (a, b) crystallographic plane. Our data do in fact contain 
relevant information on the transverse dispersion of spin excitations, which can be extracted from a quantitative 
analysis, as described in the following section. 



-I- cos^ sin^ V'2 
+ sin^ ipi cos^ tp2 



B. Data analysis 

1. Model cross section 

The general problem in interpreting inelastic neutron scattering data from powder samples is an effective loss of 
information upon spherical averaging. Indeed, the quantity of interest is the dynamic structure factor SlQ^to) that 
for each channel of spin polarization is a scalar function in 4-dimensional E — Q space. In a powder sample one 
measures the spherical average of this function, a scalar function defined in 2-dimensional space: 

5'powdGr(Q,^^) = -7- / d(f> sin d0 S {[Q sin COS (f), Q ain0 sin (f>,Q cos 9], Lo) (4) 
47r J 

The transformation S^Q^u) — s- 5'powdcr('3, is not reversible. The only way the full correlation function can be 
extracted from the experiment is by assuming some parameterized model for S{q,uj) and fitting it to the measured 
Spowdci-{Q,^^)- For PbNi2V208 and SrNi2V208 we used the following analytical form for S{q,uj), as derived below in 
the theory section: 

S(Q) = P||(Q)5||(Q) + Pi(Q)5i(Q), (5) 
2S'|| (Q) = 2 cos^ Vi cos^ V2S'|| [h, k, I) 

S*!, {h + l,k,l + 1) +511 {h + 1, k, I + 3) 

5)1 {h,k + l,l + 1) {h,k + lj + 3) 

+ 2 sin^ Vji sin^ ^28' {h + I, k + IJ + 2) , (6) 
2Sjl{Q) = 2cos^ Vi cos^ tl^2S'Ah,k,l) 

+ cos^ Vi sin^ V'2 [S'x_ ih + l,k,l + 1) +5^ (/i + 1, fc, / + 3)] 
+ sin^ cos^ V'2 [S"! ih,k + l,l + 1) +S'j^ {h, k + l,l + 3)] 

+ 2 sin^ sin^ '1P2S' {h + l,k + 1,1 + 2) , (7) 

V'l = —h, (8) 
a 

V'2 = ^k. (9) 

In these formulas the argument uj has been dropped. The phases V"! and V'2 represent the 3-D structure factor of the 
spiral-shaped spin chains, and d = 0.08a is the offset of each Ni^+-ion relative to the central axis of the corresponding 
spiral chain, along the a or b axis. The polarization factors for longitudinal (polarization along the c-axis) and 
transverse (polarization in the (a, b) plane) spin excitations are defined as: 

P\\iQ)^sin\Qrz), (10) 

Pj_(Q) = l + cos2(Q7i)- (11) 

The dynamic structure factors 6*^1 and S\ for straight (as opposed to spiral-shaped) Haldane chains are written in the 
Single-Mode Approximation (SMA)ErM: 



Zv 



Sl^{Q,uj)^ ^^ ^^J ihuj-huJuiQ)), (12) 



7v 

S'AQ,^)^j^^;j^^5{nu:-hu:^{Q)). (13) 



Here v = 2. 48 J is the spin wave velocity, and Z = 1.26 (Ref. 15). Finally, the dispersion relation for weakly coupled 
chains are given by: 
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= A| +u^sin2(7r?/2) - -^ZvJi^^ cos(7rV2) [cos(7r/i) + cos(7rfc)] (1 - cos7r//2), (14) 

[hLj±{Q)f = Al+v^sin^{TTl/2) - ^ZvJi^±cos{Trl/2) [cos(7r/i) + cos(7rfc)] (1 - cos7r;/2), (15) 

where A|| and A^ are the longitudinal and transverse intrinsic Haldane gapafor non-interacting chains, respectively. 
The splitting of the triplet in our model is caused by single-ion anisotropy:t£l 

A_L = (A) - 0.57i:», (16) 
A|| (A) -I- 1.41D. (17) 

An important parameter is the mean intrinsic Haldane gap (A) = (A|| + 2A^)/3. As can be seen from Eqs. and 
[l7| , it is, to the first order, defined by J alone: (A) « 0.41J.E£l It is also useful to define the actual gaps (excitation 
energies at the 3-D AF zone-center): 

i?L.,± = Ai-2Zz;|Jia|, (18) 
= Af-2Zz;|Ji,|||. (19) 

The spherical average of Eq. I was calculated numerically using a Monte-Carlo algorithm that also eliminates the 
(5-functions in Eqs. |l2| and |l^. The parameters were then refined by a standard least-squares routine to best-fit the 
data. To accelerate the fitting process the ASD data were binned to a rectangular 0.046 A~^x 0.1 meV resolution. 
An additional benefit of this binning is that it allowed us not to worry about resolution effects. The lower 0.5 meV 
energy transfer range, that contains the elastic-incoherent and possibly phonon scattering, was excluded from the fits. 



2. Analysis of PbNi2 V2 Os data. 

For The Pb-compound the adjustable parameters were the two intrinsic gap energies Ay and Aj^, and the doublett 
3D gap i?inin,j_- The singlet 3-D gap was fixed to ii'min.it = 1.2 meV, as determined in high-field bulk measurements. cl 
The only additional parameter was an overall scale factor. The least-squares refinement yields: A^ = 4.0 ±0.25 meV, 
A|| = 3.1 ± 0.3 meV and = 2.4 ± 0.2 meV (Ji < 0), with = 2.3. The resuhing fit is shown in Fig. |. 

Substituting the obtained values values into the expression for S'(Q,a;), performing a powder average and convoluting 
the result with the spectrometer resolution function, reproduces the measured 3-axis const-i? and const-Q scans rather 
well, as shown in solid lines in Figs. ||b,^ The relatively large of the global fit to the ASD data is to be attributed 
to systematic error, primarily due to "spurious" scattering. In particular, the data are partially contaminated by 
spurious of type "fci — > ki" and "fc/ — > kf" originating from the stronger Bragg powder lines from the sample. Areas 
that are affected by these spurious processes are shown as shaded "streaks" in Fig. Another prominent streak in 
the data is to the left from, and parallel to, the two shaded areas. The origin of this feature is not known. It is 
however resolution-limited and temperature- independent, and is thus almost certainly spurious and of non-magnetic 
origin. Considering that this type of systematic error is unavoidable in powder experiments, the obtained model fit 
to the data is quite acceptable. I-. 

The refined value for £^min,± is in very good agreement with the result of high-field measurements.lj The interchain 
coupling constants are obtained from Eqs. [l8| , p^ : Ji^± = —0.18 meV and Jij| = —0.14 meV. Note that | Ji| is larger 
by a factor of 1.5-2, compared to our previous estimate {Jj_ — 0.096 ± 0.003 meV) in Ref. ||. This discrepancy should 
be partly attributed to a difference in the definition of Ji. In our previous model each spin was coupled to 4 spins 
in adjacent chains (coordination number 4). In the present model the inter-chain coordination number is 2. This 
automatically translates into a factor of 2 for |Ji|. The negative sign of Ji indicates that this coupling is actually 
ferromagnetic. Note, however, that in our model Ji couples spins offset by c/4 along the chain axis, so that the 
effective mean field coupling between interacting chains is still antiferromagnetic, as in our previous model. The 3-D 
magnetic zone-center, where transverse dispersion is a minimum, is at (1, 1,2). From the refined Ax and A|| we can 
also get the in-chain coupling constant: J « 9.0 meV. This value is in a better agreement with the high-temperature 
susceptibility estimate J w 8.2 meV in Ref. H, than our preliminary neutron result J = 9.5 meV, from the same 
reference. The anisotropy constant D can be estimated from Aj^ and A|| using Eqs. ^,0- For PbNi2V208 we get 
D = -0.45 meV. 

To convince the reader that the Ji is in fact a relevant parameter, we have performed simulations of the dynamic 
structure factors using the same values A^ — 4.0 ±0.25 meV, Ay = 3.1 ±0.3 meV, as determined in the analysis of the 
ASD data, but with the sign of Ji reversed: Ji^± = 0.18 meV and Jij| = 0.14 meV. The resulting simulated powder 
cross section is visualized in Fig. Il^a, and is clearly very different from the inelastic intensity measured experimentally. 
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Similarly, Fig. |10| b shows a simulation with Ji^± = Ji || = meV. Ji has a particular impact on the shape of the 
constant-energy powder scans. The dashed lines in Fig. |^ are simulations for uncoupled chains. The well-defined 
intensity maximum seen in the data at \Q\ ~ 1.2 A~^, is replaced with a broad monotonous feature for uncoupled 
chains. In contrast, the interacting chain model reproduces the peak rather well. 

3. Analysis of SrNi2V2 0s data 

Similar data analysis was performed for SrNi2V2 08. in this case was fixed at zero value (see discussion in 

the theory section below). The best fit to the ASD data is obtained with = 3.9 ± 0.3 meV, A|| — 2.8 ± 0.4 meV 

and i*;^™'"^ = 2.35 ± 0.3 meV with = 2.0, and is shown in Fig. ^3. These values correspond to interchain coupling 
constants Ji^± ~ —0.18 meV and Ji y = —0.15 meV, respectively, almost identical to the corresponding PbNi2V208 
values. An important consistency check is that using the gap values to estimate J gives almost the same value as for 
PbNi2V2 08, as expected: J = 8.6 meV. The easy-axis anisotropy constant in SrNi2V2 08 is larger than in PbNi2V2 08: 
D — —0.56 meV. Simulations for constant-iJ and constant-Q scans measured for SrNi2V2 08 in the standard 3-axis 
mode are shown in solid lines in Figs. |^a,^. 

4- Observation of actual spin gaps at the 3D magnetic zone-center 

From the point of view of spin dynamics, the main distinction between the singlet-ground-state PbNi2V2 08 and the 
magnetically ordered SrNi2V2 08 is the presence of a spin gap in the former system, and a gapless excitation spectrum 
in the latter. The gap is directly accessible experimentally at the 3-D AF zone-center, where the dispersion of magnetic 
excitations is a global minimum. As mentioned above, for both vanadates the 3D zone-center is at Q^"-* = (1, 1,2), 
which corresponds to a momentum transfer IQ^"^] ~ 1.67 A~^. Constant-Q scans extracted from our ASD data for 
this momentum transfer, as well as a standard Ef = 5 meV constant-Q scan at this wave vector, are shown in Fig. 
To better understand these data, we note that near the 3D zone-center the single-mode contribution to the dynamic 
structure factor, independently of the details of the spin Hamiltonian, should be of the following form: 

S{Q,uj)(x—5{lu-ujq), (20) 

[huQf = + «^(Q|| - q[|"V + ^l{Qi- QTf- 

In this formula A is the gap energy, an d vh and v± are spin wave velocities along and perpendicular to the chain axis 
c, respectively. The data shown in Fig. |ll| c were collected with a horizontally-focusing analyzer. At energy transfers 
below 1 meV in this mode we are picking up a great deal of diffuse and phonon scattering ( (1, 1,2) is an allowed 
nuclear Bragg peak !). Above this contaminated region though, in a powder sample we effectively observe a Q-integral 
of the cross-section around the 3D zone-center. The same applies to the scans in Figs. pT| a and b, where integration 
was performed in the range 1.6-1.7 A^^. For the Q-integrated intensity Eq. |2^ gives: 

Siuj) cx Vi^^ - A2. (21) 

The linear increase of intensity seen for SrNi2V2 08 in Fig. ^ should thus indeed be interpreted as due to a gapless 
spin wave (A = 0). In contrast, for PbNi2V208, the threshold behavior is a clear sign of an energy gap. These effects 
are also seen in constant-Q scans extracted from the ASD data by binning the pixels in a 0.1 A-i U-range. 

C. SrNi2V208: Temperature dependence 

To better understand the mechanism of long-range ordering of SrNi2V208 we studied the temperature dependence 
of inelastic scattering at the momentum transfer = 1.67 A~^ in this compound. Typical data are shown in 

Fig. [l|. The solid hues were obtained by fitting our model cross section to the data at each temperature. The lower 
0.75 meV energy transfer range was excluded from this analysis, as the difficult-to-estimate phonon contribution in 
this range is expected to increase dramatically with increasing T. The set of independent parameters was slightly 
modified. The anisotropy splitting of the triplet D, as well as both inter-chain coupling constants Ji ^ and || were 
fixed at the values determined at T = 1.5 K. The mean gap (A) and an intensity prefactor were refined to best-fit 
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the scans at each temperature. The obtained temperature dependences are shown in Fig. |T^. As previously observed 
in other Haldane-gap systems, (A) increases with increasing T. Even though this change is rather small, according 
to Eqs. ^,|9[ it corresponds to an appreciable variation of the gap in the longitudinal mode (Fig. Upon cooling, 
the longitudinal gap approaches zero at T = Tn, which results in a |SQ£t-mode transition to a magnetically ordered 
state. This type of behavior is very similar to that found in CsNiCl3.laQ 



V. THEORY AND DISCUSSION 



A. Derivation of the model cross section 



To calculate the effect of weak inter-chain interactions-pn dynamic spin correlations in PbNi2V208 and SrNi2V208 
we shall use the Random Phase Approximation (RPA) .llZl This approach for directly coupled Haldane spin chains has 
been successfully applied, for exanaple, to CsNiCls (Ref. ^). More recently it was shown to also work well for Haldane 
chains coupled via classical spinsJiS The only additional difficulties in the present case arise from the complicated 3-D 
arrangement of magnetic ions in the PbNi2V208 structure, and from the rather non-trivial geometry of interchain 
coupling. To somewhat simplify the task we shall break it up into two distinct problems. First, we shall worry only 
about the topology of magnetic interactions and consider an equivalent Bravais lattice of spins, assuming straight spin 
chains and a system of inter-chain bonds shown in Fig. ^jb. Making use of the general RPA equations for coupled spin 
chains, summarized in Appendix I, we shall write down the RPA susceptibility and dynamic structure factor S'{Q, lo) 
for this model. In a separate step we shall adapt the result to the more complex structure of PbNi2V208, using the 
formulas of Appendix II. 



1. Structure factor for an equivalent Bravais lattice 



The first step in the RPA calculation is to write down the bare (non-interacting) susceptibilito- foj- an isolated 
Haldane spin chain, for which the single-mode approximation (SMA) is known to work rather well:E3iij 

1 - cos(g3c/4) Zv 

Xo[Q,Uj) = .2 I 2 ■ 2.^ 77 -—^ (22) 

This expression should be used separately for each channel of spin polarization, with appropriate values of gap energy 
A for each particular mode. 

Expression is to be substituted in the general RPA equations Even for the "straightened out" spin chains in 
Fig. ^ there are still 16 spins per unit cell, so we end up with 16 couple equations and 16 modes for each polarization! 
For the low-energy part of the excitation spectrum it is however quite appropriate to use the approximation of Eqs. [sj. 
This reduces the problem to only 4 self-consistent RPA equations (there are 4 spin chains in each crystallographic 
unit cell). Moreover, at this level of approximation the set of 4 RPA equations is degenerate, and we end up with a 
single RPA equation for each polarization: 

Xrpa(Q) - Xo\Q) [xo{Q)J{Q) + 1] , (23) 

where 



J{Q) = Ji cos(Q.c/4) [cos{Q^a/2) + cos{Qya/2)] . (24) 

This equation is easily solved analytically. By taking the imaginary part of the thus obtained Xrpa(Q) we arrive at 
Eqs. |l| |l|. 



2. Actual structure 



The transformation from the Bravais spin lattice (straight chains) to the spiral-chain structure of PbNi2V2 08 and 
SrNi2V2 08 is schematized in Fig. |l5|. This figure also serves as an illustration of the notation used in Appendix I. 
Applying Eq. H4 is straightforward, and we will spare the reader the tedious calculations. The final relation between 
S{q,LL!) and 5*^,0;) is as in Eqs. 
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B. Temperature effects 



1. Disordered phase 

The derivation above can be repeated almost verbatim for the case T > 0. Temperature enters RPA calculation 
indirectly, through an intrinsic temperature dependence of bare susceptibilities of individual chains (Eq. ^). Two 
effects at T > need be considered: i) the increase of the Haldane gap, compared to its T = value, and ii) damping of 
Haldane excitations. Both phenomena have been observed in a number of model quasi-l-D systems (see, for example, 
Refs. p^-pl|). Unfortunately, it is almost impossible to extract meaningful information regarding excitation lifetimes 
from powder inelastic data. As most of the observed inelastic signal originates from excitations with energies greater 
than 2 meV, the temperature dependence of the intrinsic Haldane gap is much easier to observe. For this reason 
we analyzed the temperature dependence of the inelastic signal measured in SrNi2V208 using the same single-mode 
approximation (E g. p2| ) , as at base temperature (see fits in Fig. |l2|) , to obtain the T-dependences of intrinsic Haldane 
gap energies (Fig. |13D. From the mapping of the Heisenberg Hamiltonian on the quantum non- linear a- model (NLSM) 
one expects a rather steep increase of the gap energy with increasing temperature:E3 

A(0). 



A(r) « A(0) + V2^yTA(0) exp(--^). (25) 

The actual increase of the gap energy observed experimentally NENP, NINAZ,0 and YaBaNiOg,!] was found to 
be consistently smaller that this NLSM prediction. The same discrepancy, namely a rather slow increase of the gap 



energy with temperature, is seen in SrNi2V2 08 as well. Equation is plotted in a dashed line in Fig. 13 a,. The 
solid line is an empirical fit in which the prefactor \/27r in Eq. ^ was replaced by an adjustable parameter A. For 
SrNi2V208 we get A — 0.67(1) ^ ^/2n. A self-consistency check for our SMA model is that the refined intensity 
prefactor is almost T — independent, as shown in Fig. [l3|b. In other words, the decrease of actual excitation intensity 
is entirely due to the increase of gap energy, and the intensity scales as l/uj. 

In SrNi2V2 08 the ordering temperature Tn = 7 K is significantly smaller than the intrtinsic Haldane gap energy 
A « 4 meV. As both the gap energy and excitation width are expected to increase exponentially with T, such a 
small T/v suggests that inter-chain interactions in SrNi2V208 are barely strong enough to produce a LRO-ground 
state. Long-range ordering in SrNi2V208, and the absence of such in PbNi2V208, are a very "lucky" coincidence that 
results from a fine interplay between in-chain interaction strength, magnetic anisotropy and inter-chain coupling. 



S. Magnetically ordered state 

Our model cross section was derived under the implicit assumption that the system is in a non-magnetic state, 
which is not applicable to SrNi2V2 08 at T < Tn. In this regime one expects an increase of the gap energies, at least 
for the longitudinal rnpde- due to the presence of a mean static staggered field generated by the ordered staggered 
moment in the system.c2lc3 To adapt the MF-RPA calculation to this temperature regime one has to know the ordered 
moment in the system. To date, powder diffractiop-|experiments on SrNi2V208 failed to detect any magnetic Bragg 
rcfiections in this compound at low temperatures,E3 and the magnetic order parameter is thus expected to be very 
small. This is consistent with our previous observation that the system is "barely" 3-dimensional enough to become 
ordered at low temperatures. Provided the ordered moment is small, within the accuracy of our powder experiments 
it seems quite appropriate to use the same cross section for T < Tn, as for the paramagnetic phase, and postulate 
-^min,|| = 0- This assumption is equivalent to assuming Tn — 0. 



3. Placement on the phase 

The line of quantum phase transition separating the ordered and spin-liquid states can be derived from Eqs. |l6-19. 
The critical value for Ji corresponds to the lower gap energy (in our case i?min,||) being equal to zero. The resulting 
phase diagram is very similar to the direct numerical calculations by Sakai and Takashi,!] and is shown in Fig. |l| Using 
the results for PbNi2V208 and SrNi2V2 08 described above, we are able to place the two new materials on the same 
plot. It has to be emphasized, that for SrNi2V208 the parameters D and Ji were extracted from the experimental 
data using a cross-section for a disordered system, and assuming -Emin.ll = 0. That SrNi2V208 lands exactly on the 
phase boundary is thus an artifact of our data analysis procedure. In reality, SrNi2V2 08 must be positioned slightly 
above the phase boundary. The disordered PbNi2V208 is not infiuenced by such an artifact, and its positioning can 
be considered quite reliable. 
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VI. SUMMARY 



Our neutron scattering results help explain the different ground state properties of the two very similar vanadates. 
The ratio of inter-chain to in-chain coupling in SrNi2V208 is slightly larger than in PbNi2V208. SrNi2V208 is also 
driven towards LRO by the somewhat larger magnetic anisotropy. The quantitative data analysis enables us to 
precisely position the two compounds on the Sakai-Takahashi phase diagram, just opposite to each other relative to 
the ordered-disordered phase boundary. This almost unbelievable coincidence opens many exciting possibilities for 
future studies. Experiments on aligned powders will enable more accurate measuijements of the excitation spectrum 
near the 3D zone-center. High-pressure studies, as previously attempted for NENP,E3 may lead to the first observation 
of pressure-induced long-range ordering in a quantum-disordered magnet. 
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APPENDIX I: RPA SUSCEPTIBILITY FOR WEAKLY COUPLED SPIN CHAINS 

The general RPA technique of calculating the dynamic susceptibility and structure factor for weakly interacting 
systems is well established and documented (see, for examj}le_Ref. ^). In particular, the approach has been on many 
occasions applied to weakly coupled quantum spin chains O'Q In the present section, without any claims of novelty, 
and for reference only, we shall derive some useful RPA results for weakly coupled spin chains assuming a rather 
general geometry of inter-chain bonding. Note that the notation used in this Appendix is totally independent from 
that in the rest of the paper. 

We consider a crystal structure composed of identical parallel uniform spin chains that run along the a axis. The 
spacing between spins in each chain is a/K, where K is integer. There are M chains in each crystallographic unit 
cell. The origin of the m-th chain is at ii„ -I- pm- The position of the i-th spin in chain (n,m) can be written as 
Rn + pm + ia/K. A single crystallographic unit cell thus contains MK spins. For convenience we shall break up 
the in-chain spin index i into two: i = IK + k, < k < K . The spin Hamiltonian for inter-chain interactions will be 
written as: 

iJinterchain = E E E E E ^ll' '"^ 4' i"'^ -^(^^ ^' ^ ^1 + fc, Kl' + k') (26) 

n,n' m.m' 1,1' k—0 k' —0 

We assume all exchange interactions to be simultaneously diagonal in spin projection indexes. All the equations in 
this section thus apply to a particular channel of spin polarization. In our model the exchange constant satisfies 
certain translational-symmetry relations: 

J(n, n', TO, to', Kl + k, Kl' + k') = J{{Rn' - Rn), m, m', (l' - /), k, k'). (27) 

In order to write down self-consistent RPA equations, we have to answer the following question. Suppose we have 
artificially induced a spin density s^™ \r) (or, equivalently, s^™ •'(q'))in the chain-sublattice to'. What exchange field 

will this spin density project on the chains in sublattice to? By definition, the exchange field acting on spin s|'^'™'' is 
given by: 

M-l K-1 

/^lr^ = E E EE4"l'"'^^((«n'-^n),-,-'a'-o,fc,fc') (28) 

n' ni'—O I' k'—O 

We, of course, will be interested in the Fourier transform of the exchange field acting on sublattice to: 
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K-1 



h^^\q) [-^^(Rn + P.n + a{k/K + 0)] . 



n I k=0 



Having introduced the definition 

Jk,m,m'{q) = ^ cxp [-iq{Rn + p,n' ~ Pm + al + a{k' - k)/K)] J{Rn,m',m,l, k' , k), 



(29) 



(30) 



.,k'.l 



it is straightforward to verify that: 



A/-1 K-l 



= K=0 



where 



JK,7n,m'{q) = ^ X! Jk,m,m'{q) exp [Zmnk / K] . 



(31) 



(32) 



Equation ^ enables us to immediately write down the self-consistent RPA equations: 

-1 



M-l K-1 



X! 51 ^o" (9 + «"*'^)'7K,m,m'(9) + 1 



_m'— K— 



(33) 



Here Xo^^\qi^) f^^id x'-™'(q,a;) are the bare (non-interacting) and RPA-corrected wave vector dependent dynamic 
susceptibilities of spin chains in the rn-th sublattice, respectively. For any wave vector q one obtains KM equations: 
susceptibilities at K wave vectors become coupled to each other and there are M independent chain-sublattices. 

For a Heisenberg AF chain, at energy transfers much smaller than the in-chain exchange constant, Xo™'(9) 
typically very small except when momentum transfer along the chain axis is close to qo, defined by qoa = tt. To a 
good approximation we can thus replace Eq. B3^ with: 



X^"^\qa^7r,u;) 
X^^\qa^7r,u;)^0 



M 



-, -1 



\'li^)Jm,m'iq) + 1 



.m' — 1 



(34) 



where 



K K 



Jm,m'{q) = J^YYYY J{Rn,m,m',l,k, fc')exp [-iq(Rn + al + a{k' - k)/K + {p,n' - Pm))] (35) 

n I k=lk' = l 

is simply the Fourier transform of exchange interactions between sublattices m and m' . Note that there are only AI 



equations in the system 34 



APPENDIX II: EQUIVALENT BRAVAIS LATTICE 



In this section we shall derive a useful relation between magnetic dynamic structure factors of two spin lattices 
described by the same Hamiltonian, but featuring different 3D arrangements of magnetic sites. Lattice I is assumed 
to be a simple Bravais lattice of magnetic sites and consists of N unit cells of volume v. The origin of the n-th unit 
cell Rn'' coincides with the position of the n-th spin: ri^-* 



Rn ^ ■ Unlike lattice I, lattice H is non-Bravais, having M 



spins per unit cell, and a unit cell volume V—Mv. There are K = N/M unit cells with origins at R^]^^\ respectively. 
Lattice H is obtained from lattice I by shifting each spin n by p„. In lattice II the spins are thus positioned at 



Pn- The relation between the two lattices allows us to write r„ 



(11) 



R 



(11) 



Jii) 



R 



(II) 



r'r. 

( n = kM + m, < m < M). We assume the spin systems to be equivalent in the sense that the spin Hamiltonian 
of system II written in terms of site-spin operators s„ is identical to that of system I. The notation introduced above 
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is independent of that used in the rest of the paper or in Appendix I, and is illustrated, for the particular case of 
PbNiaVzOg, in Fig. p. 

The Fourier transform of the total spin operator for lattice I is defined as: 



N-l K-lM-1 

s^'Hq) = ^ s„ exp[-iqi2(^)] '^«+rn exphiq(i2i5, + r(„^))]. (36) 

Again, for simplicity, we have omitted the spin projection indexes. The reverse relation is given by: 

SkM+,n = J dqs^'\q) exp[iq{Ri'l + r^^))], (37) 

where the integral is taken over the Brillouin zone of the Bravais lattice ( a reciprocal-space volume K times as large 
as the Brillouin zone of the lattice II). Similarly, for lattice II, by definition: 

K-l M-l 

s^''\q) = E E ^"'^+™ exp[-iq(i2f ) + r^^^))]. (38) 

fc=0 m=0 

The equivalence of the two systems allows us to directly combine the last two equations: 

K-l M-l 

^^''^9) = E E eM-iq{Rl'"+rlrV^)] I dpeMHRl'li+rliW^ip). (39) 

k=0 m=0 



9) - 7^ E E exp[-iq(fii^^' +r(f ))] / dpeM^p{Ri2, + rg^)]s('\p). 
The expression can be simplified by noting that 

E exphi(g - p)Ri"^] =^-^Y.^[p-q-T). (40) 



V 

k=0 



Here the inner sum is taken over reciprocal-lattice points T^, k ~ 0...K, for lattice II. Plugging this result into the 
expression into Eq. gives: 



M-l 



s^''Hq) = ^ fdP ^^'^(P) E -^(P - 9 - ^) E exp[-i(gr-(„") - prii^)] 

T m=l 

= M E exp(-i9P™) E ^''^(9 + ^m) exp(iT^r(^)) (41) 

m— M— 1 

In this the last expression the first sum is taken over all reciprocal-lattice vectors of lattice II within the first 
Brillouin zone of lattice I. 

Equation ^ for the Fourier transform of spin operators can now be directly plugged into the definition of dynamic 
structure factors for the two lattices: 

{q, dt exp(-ic.O (s(^) (-q, O)s(^) (q, t))) (42) 

S^"\q,uj) I diexp(-ic.t)(s(^^)(-q,0)s(")(q,i))). (43) 

2ttTi J 

A straightforward, though somewhat tedious calculation gives the result that we were after: 

2 



M-l 



5(^^)(q,c.)=E 
11=0 



M-l 



E exp(-iqp,„)cxp(iT^r(^)) 



m=0 



S^'\q + T^,uj) (44) 
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FIG. 1. D — Ji Phase diagram for weakly coupled Haldane spin chains,u showing the location of some well-characterized 
quasi-lD S — 1 systems. For each material Ji is scaled by the coordination number. CsNiCla data from Ref. ^ Y2BaNi05 
data from Ref. H; NENP data from Ref. ||; AgVP2S6 data from Ref. PbNi2V208 and SrNi2V208 are placed based on the 
results of this work. 



FIG. 2. (a) Crystal structure on PbNi2V208 and SrNi2V208. The magnetic chains are formed by edge-sharing NiOe 
octahedra (dark gray) and are bridged by VO4 tetrahedra (light gray), (b) A perspective view of the spiral-shaped Ni^^-chains 
in PbNi2V208 and SrNi2V208. 



FIG. 3. Geometry of inter-chain interactions in PbNi2V208 and SrNi2V208 (dark bonds). The strongest inter-chain coupling 
is assumed to be between nearest-neighbor Ni^+ ions from adjacent chains. These sites are offset relative to each other by c/4 
along the chain axis. 



FIG. 4. Topology of magnetic interactions PbNi2V2C)s and SrNi2V208. Inter-chain bonds are shown as dashed lines. Solid 
lines are in-chain Ni-Ni bonds, (a): Model used in Ref. H. (b): Model used in the present paper. 



FIG. 5. False-color plot of inelastic neutron scattering intensities measured in PbNi2V208 (a) and SrNi2V208 (b) at 
T = 1.5 K using an area-sensitive detector setup. The background has been subtracted, as described in the text. Shaded 
areas indicate the location of possible spurious. 
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FIG. 6. Constant-energy scans measured in PbNi2V208 (a) and SrNi2V208 (b) powder samples at T = 2 K in the standard 
3-axis mode. The solid lines are simulations based on parameters obtained in a global fit to the ASD data, as described in the 
text. Dashed lines are similar simulations for non-interacting spiral-shaped spin chains. 

FIG. 7. Typical constant-Q scans measured in PbNi2V20g powder sample at T = 2 K in the standard 3-axis mode. The 
solid lines are as in Fig. ^. 

FIG. 8. Typical constant-Q scans measured in SrNi2V208 powder sample at T = 2 K in the standard 3-axis mode. The 
solid lines are as in Fig. ^. 

FIG. 9. Simulated inelastic scattering cross section based on parameters obtained in a global fit to the ASD data for 
PbNi2V208 (a) and SrNi2V208 (b). 

FIG. 10. Simulated inelastic scattering cross section based on parameters obtained in a global fit to the ASD data for 
PbNi2V208, assuming a reversed sign of Ji (a) and no inter-chain interactions (b). 

FIG. 11. Constant-Q scans extracted from the ASD data sets for PbNi2V208 (a) and SrNi2V208 (b) at a momentum transfer 
that corresponds to the 3D antiferromagnetic zone-center (symbols). Unlike in SrNi2V208, the measured inelastic intensity in 
PbNi2V208 extrapolates to zero at a non-zero energy transfer (arrows). This is even better seen in constant-Q scans measured 
in the standard 3-axis mode (c). Open and solid symbols correspond to SrNi2V208 and PbNi2V208, respectively. In the latter 
plot the data for the two compounds were brought to a single scale using the elastic incoherent signal. In all cases the solid 
lines are as in Fig. ^ 

FIG. 12. Typical constant-Q scans measured in SrNi2V208 powder at different temperatures, at a momentum transfer that 
corresponds to the 3D antiferromagnetic zone-center. The solid lines are fits to the model cross section, as described in the 
text. 

FIG. 13. Measured temperature dependence of the mean intrinsic Haldane gap (A) (a) and intensity prefactor in the model 
cross section described in the text (b) for SrNi2V208. 



FIG. 14. Measured temperature dependence of the 3D gap energies in for SrNi2V208. The solid lines are as in Fig. 

FIG. 15. Relation between the actual structure in PbNi2V208 (b) and the equivalent Bravais lattice of spins (a). The 
notation is as in Appendix II. 
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